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Abstract 

Symmetries affected by the anomaly do not survive quantization and cannot be un- 
derstood classically. They are of fundamental importance and offer an opportunity of 
expanding the theoretical framework. We examine the theory of the anomalous sector, 
starting from lowest order Chiral Perturbation Theory (ChPT), leading up to the con- 
struction of the recently developed Lagrangian of O (p^) describing anomalous processes. 
This Lagrangian contains a set of chiral coefficients that must be determined phenomeno- 
logically. Using currently available experimental data, we fit as many of these coefficients 
as possible. The results of the ChPT treatment are then used to test the validity of the 
two main alternative models employed in the anomalous sector - Vector Meson Dominance 
(VMD) and the constituent Chiral Quark Model (CQM). 
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1 Introduction 

In high energy QCD, we have asymptotic freedom and scattering can be calculated in a power 
series expansion in the strong coupling constant. The non-perturbative effects are determined 
phenomenologically in the form of structure and fragmentation functions. Predictions take 
the form of relations between amplitudes parametrized by ag and the structure functions. 

In low energy QCD, the predictions are relations between amplitudes with a structure 
determined by symmetry constraints, parametrized by empirically determined coefficients. 
The power series expansion in is replaced by an expansion in the low energy. Heavy fields 
of the high energy region can be integrated out and their effect is consequently encoded in 
the coefficients of the low energy effective theory. Proceeding up to O (p^) in ChPT, we 
are confronted with a large number of empirical parameters. In principle, QCD should be 
able to predict these parameters, as well as the structure functions appearing in the chiral 
Lagrangian. No rigorous derivation exists up to date. The combined use of tentative models 
and phenomenological knowledge provides us with insight of the physics leading to the chiral 
Lagrangian. The aim of this thesis is to solve for as many low energy chiral coefficients as 
possible by making use of the available experimental data in the form of widths, slopes and 
form factors. 

Classically conserved currents that are affected by the anomaly do not survive quanti- 
zation. To evaluate the ensuing effects, it is necessary to employ a more rigorous analysis, 
taking quantum corrections into account. The effect of the anomaly in the chiral Lagrangian 
framework was first analyzed by Wess and Zumino who realized that the result could 
not be expressed in a local effective Lagrangian. Their result, in the form of a Taylor series 
expansion, was later given an elegant geometrical interpretation by Witten 0. 

So far most articles treating anomalous processes quote the leading order amplitudes and 
meson one-loop corrections. Higher order contributions are then estimated using models like 
Vector Meson Dominance (VMD) or the more recently developed chiral Constituent Quark 
Model (CQM). In this thesis we will instead test the validity of VMD and CQM by comparing 
with the results of the experimentally fixed chiral O (p^) coefficients. 
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2 Symmetries 



2.1 Chiral Symmetry 

Apart from the more obvious symmetries of the standard model, such as the gauge symmetries 
SU (3)^ X SU {2)j^ X U {1)y , there exist global vector symmetries like the fermion number 
symmetries, the isospin symmetry, and the less accurate SU (3) flavor symmetry. The flavor 
symmetries are valid if the masses of the quarks included by the symmetry are set equal. 
Therefore, the isospin symmetry is broken by the up-down quark mass difference, as well as 
by electromagnetic and weak interactions. By imposing the considerably stricter condition, 
niq = 0, so called chiral symmetries arise. Since there are no longer any mass terms to 
couple fields of differing chirality, the fields are invariant under separate left- and right-handed 
transformations. The QCD Lagrangian in the massless limit is given by 

^Qd% = E Ql" (id, - gs^Gl) q - IG^^^G^'' 

q=u,d,s 

= i,Lpi^L + i^Rp^R-\G%G>^'' (1) 

where V'-R.-L = ^ (1 =t 75) ip = Pr^l^P are the chiral projections of the column vector containing 
the relevant quark fields {u,d for chiral SU (2) and u,d,s for chiral SU (3)). Taking m„ = 
nid = 0, the QCD Lagrangian is invariant under chiral rotations of the fields 

V'L ^ e-'^^-^i^L = gii^L (2) 
^ e-''''-^iPR = gRijR (3) 

where {r* = ^o"*} {i = 1,2,3) are the Pauli matrices and {O^r} are the components of an 
arbitrary constant vector. The chiral SU (2) invariance is referred to as the direct product of 
the left- and right-handed chiral transformation groups, SU (2)^ x SU (2)^ . This symmetry 
is easily extended to include the strange quark, substituting the Pauli matrices for the Gell- 
Mann 3x3 SU (3) matrices |{A"} (a = 1,...,8). In eq. (0), we also find the vector U {1)y 
and the axial U invariances 



ip (4) 
^ ^ e^^-z^o^ (5) 

So apart from the normal SU (3)^ color gauge and discrete symmetries, the chiral QCD 
Lagrangian also has the global invariance 

U{l)^xU{l)yxSU{3)j,xSU{3)L (6) 

where the there are now separate invariances under SU (3)^^ and SU (3)^ , for the three 
massless quarks. The U (1)^ is conserved, ensuring baryon number conservation. The axial 
[/ (1)^ current is not conserved in the full quantum theory, and is explicitly broken by the 
Abelian anomaly. 

Let G be the direct product group of the left- and right-handed chiral transformations: 
'>pL,R — *■ gL,R'4'L,R, 9l,r £ G. Then ^ is invariant under G. However, the QCD vacuum 
structure does not share this invariance. The effect of a non-vanishing quark condensate. 
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i.e. a non-zero vacuum expectation value (0| ipip |0) 7^ 0, is the spontaneous breaking of the 
axial part of G. Spontaneous (dynamic) symmetry breaking is said to occur whenever the 
symmetry of the Lagrangian is not shared by the ground (vacuum) state. For G we have 

SU (3)^ X SU (3)^ ^ SU {3)y^L+B. (7) 

where SU (3)y = H is the remaining unbroken subgroup of simultaneous transformations 
of differing chirality. According to the Goldstone theorem, whenever a global continuous 
symmetry is broken, massless so called Goldstone (GS) bosons appear. Since it is the axial 
part that is broken, the resulting GS bosons are pseudoscalars. An SU (n) matrix has n? — 1 
parameters (due to the unitarity and determinant constraints). The matrix of the coset 
space G/H = SU {3)j^ x SU {>i)j^ / SU (3)y , used to parametrize the GS bosons is an SU (3) 
isomorphism, also represented by a 3 x 3 matrix, containing 8 pseudoscalars. These are 
the light pseudoscalar mesons of the low energy spectrum: ,7r^ , , , and t]. If we 
instead would have chosen to focus on the up and down quark, setting m„ = = 0, the 
resulting spectrum would have been described by the 2x2 SU (2) matrices, containing the 
2^ — 1 = 3 pions: tt^,tt~ and vr'^. 

The only visible symmetry in hadronic states is U (l)y x SU {3)y . The effect of U (l)y 
can only be seen by including the (non-trivially transforming) baryon fields in the effective 
Lagrangian. 

The proof of the breaking of the axial symmetry lies in the predictions it makes, and can 
be produced by Monte Carlo lattice gauge techniques. The main motivation for the existence 
of this mechanism lies in its phenomenological success and theoretical consistency. 

GS bosons are massless, but the physical particles are clearly not. However, the quark 
masses are small compared to the breaking scale of the chiral symmetry {mu^d,s < ~ 
IGeV). This enables us to treat them as a small perturbation, and we can include the quark 
masses in the Lagrangian as external scalar fields. The inclusion of mass terms will mean 
the explicit breaking of G. Perturbational treatment is supported by the relatively small size 
of the pseudoscalar masses compared to the hadronic scale (m^/m^ ~ 0.03 for SU {2) and 
mj^/ml ~ 0.3 for SU (3)). 

The task at hand is to construct a chiral invariant low energy effective quantum theory 
with the GS bosons as dynamic fields. 

2.2 Non-linear Realization of Chiral Symmetry 

Before we can construct an invariant Lagrangian, we must examine the transformation prop- 
erties of the Goldstone boson degrees of freedom. G is a compact connected Lie group, that 
is dynamically broken into the subgroup H. The coordinates of the coset space G/H = 
SU (3)^ X SU {3)j^/SU (3)^_(_^ are the remaining degrees of freedom describing the system, 
i.e. the GS bosons. The transformation of an element parametrized by the GS bosons under 
G 

gu{i^)^u{i[)h{iug) ;geG,heH 

is non- linear in nature since the generators of H and G/H do not commute. What we want 
is to construct operators that transform linearly under G, that is to say are chiral invariant. 
To this we then add the explicitly symmetry breaking light quark masses. Linear operators 
can be constructed by starting from projections onto the transformation subgroup H. As 
shown by Callan, Coleman, Wess & Zumino _4', this is the most general way of constructing 
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operators linear under G in terms of the GS bosons. A given field tp, transforming linearly 
under H, transforms under G as — > h{g,(^i)7ph~^ (9,^,1) ■ This means that all products of 
type (u, u^^ ■ if; ■ {u, n^) will transform linearly under G. For a more technical treatment of 
how this can be done the reader is referred to 

From the derivative of u and we can define the Hermitian operator and the covariant 
derivative V^, both of which transform linearly under G. 

= i {v)d^j,u - udfj^v)^ (8) 

V;,V = d^i^-[T^,i;] ■T^ = \(u^d^u + ud^u^) (9) 

These can be used as building blocks to construct the linear operators. To get the correct 
behavior in the low energy limit, it is necessary to have derivative interactions. This is not a 
problem since the Goldstone fields can easily be rewritten to accommodate this soft behavior. 
The invariant term is realized by taking the trace in flavor space, denoted by (...) . 

The parametrization of the coset space is not unique. It can be shown that S-matrix 
elements are invariant under a field redefinition -0 = (pF {(j)) , where -F (0) = 1. The power 
counting depends on the number of derivatives, which is unchanged by the redefinition, so 
the results are representation independent, order by order. It is advantageous to use an 
exponential parametrization in 3 x 3 flavor space. 

u' = U = e^^*^ (10) 
where F is a dimensional constant of same dimension as M - the Goldstone boson matrix 



i=l 



where Aj are the Gell-Mann SU (3) matrices, (j)^ the Goldstone fields and rj^ denotes the octet 
component of ry (see section on r/g-jyo mixing). 



3 Effective Quantum Field Theories (EQFT) 

QCD is a gauge theory describing the interactions of quarks and gluons. Below it becomes 
highly non-perturbative. As a consequence, we can no longer use the partonic degrees of free- 
dom to describe it. The effective chiral Lagrangian approach makes use of the very simple low 
energy spectrum of light pseudoscalar mesons: -k^ ,tt^ , , , and 77, as dynamic fields 
describing the theory. Considering the weak nature of interactions amongst and between the 
light mesons and nucleons, the perturbative approach can be reinstated by simply transform- 
ing to the relevant degrees of freedom, i.e. the light pseudoscalar mesons. This is the principle 
for all effective quantum field theories. At a given energy, only certain degrees of freedom 
are relevant in describing the theory. The non-relevant degrees of freedom can be integrated 
out, and their effect is consequently encoded in the coefficients of the effective theory. All 
quantum field theories can be regarded as EQFTs. However, they differ in their degree of 
renormalizability. Since the theory is valid only below a given intrinsic scale parameter A, 
we can expand amplitudes in terms of E/ and require a finite number of counterterms in 
order to regularize at any O (E^/A"') . 
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If we want to be able to study relevant physical processes, the needed external source 
fields must enter the Lagrangian. These are conveniently included in a chiral invariant way. 
They take the form of external scalar (s), pseudoscalar (p), right- (r^) and left-handed (£^) 
3x3 matrix source functions. 

^QCD = ... - ip^^iPL^^ip - tp-y^iPRr^^J -i'L{s + ip) ipR -ipRis- ip) IpL (12) 

Here we see that the quark mass matrix s — tti — d%cig{TTiui md,ms) exphcitly breaks the left 
and right chiral symmetries. The electroweak gauge fields enter through and r^. 

The low energy effective action for the GS bosons is a functional of external sources. We 
obtain the QCD connection by considering the effect of the sources. 



Only the GS bosons are the relevant degrees of freedom at low energies, so in integrating out 
the heavy fields (thereby absorbing them into coefficients) and transforming to the GS system 
gives 



4 Chiral Perturbation Theory (ChPT) 

Chiral Perturbation Theory is an EQFT describing hadronic interactions in the low energy 
limit of the standard model. It is valid below the breaking scale of chiral symmetry, i.e. for 
energies <C IGeV ~ A^. The chiral theory successfully describes the mesonic sector, and has 
also been extended to other fields such as heavy quark and bound state dynamics. ChPT is 
the evolved form of Partial Conservation of the Axial current (PCAC) and current algebra 
techniques. 

4.1 Lowest Order Effective Lagrangian 

To obtain terms that are invariant under both chiral and Lorenz symmetries, at least two 
derivatives of u or must be present. In the absence of external fields there is only one term 
of O 

= (d^Ud^'uA (15) 



-2 

where the coupling is set to reproduce the correct kinetic term. Adding source functions 

coupled to their associated currents (following Gasser & Leutwyler [S]), as in eq. ((T^ . we find 
the global symmetries © implying the following invariances: 

V' — > e'^ip (16) 

lb — > e'^^^i) (17) 

= V/^-af,—^ giif^al^ (19) 

r,^ = V/i + a/i—^ 9Rri,gR^ (20) 

s±ip gR^L{s±ip)gl\ (21) 
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These can all be made local by adding terms to the transformation of the vector fields 



^ i 9L {l^ + id^)gl^ 

The pseudoscalar source is not relevant for the processes considered here (it appears for 
example in the Higgs-sector) , and the scalar field is set to the light quark mass matrix. 
(12) shows that the s field transforms non-trivially under G, explicitly breaking the chiral 
symmetry. But because the light quark masses are much smaller than the chiral breaking 
scale, we still have approximate chiral symmetry. 

Comparing with the SM Lagrangian we can make the identifications 

v^. -eQA^-^{T+W; + T^W^) (23) 

^ 'T+W+ + T-W~) (24) 



2V2 




where the neutral weak gauge fields have been omitted, since their contribution will be strongly 
suppressed by the heavy Z*^ mass. Q is the electric charge matrix Q = \diag (2, —1, —1) and 



(25) 

and its Hermitian conjugate T_ contain elements from the weak mixing matrix. 
The non-abelian field strength tensors are given by 

p^^ = ^^'t -d'^t -i[e,t] (26) 
F^^ = d^'r" -d^r^" -ilr^'^r"] (27) 

Local invariance is maintained by replacing in (|15jl by the covariant derivative D^. 

d^U ^ D^U = d^U - ir^U + iUl^ (28) 



4.2 Power Counting 

To organize our results we must examine how powers in the chiral expansion are counted, so 
that we can assign them to the Lagrangian of appropriate order. It is the order expansion 
that makes the whole approach useful. If we assign d^u and a^, the same power counting, 
D^U becomes a first order homogeneous term in the derivative expansion. We then have 
U ^ O {p^) , afj_,v^, ~ O (p^) , s,p, -F^^ ~ O (p^) . The lowest order Lagrangian including 
external sources is given by: 

J^2 = ^ [d^UD^U^ + xU^ + Ux^) (29) 

where x = 2i?o (s + ip), and Bq is a constant related to the vacuum expectation value (^ijjip)^. 
The coupling F is can be identified with the pion decay constant Ft^. 
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4.3 Vacuum Expectation Values and Masses 

The axial symmetry is a hidden symmetry, meaning that it is dynamically broken - the 
invariance of the Lagrangian is not shared by the vacuum state. It can be seen that the 
vacuum state does not transform separately under left and right chiral transformations, as it 
couples ^pL with tpR. 

(oiviV'lo) = {o\i/j{PL + PRfi^\o) 

= {0\i^L^R\0) + {0\ijR^L\0) (30) 

where we have used the projection operator properties and 

4^L,R = {PL,Ri^)^ 7° = V'^^L.ijT^ = i^PR,L (31) 

(by the commutation relations for the gamma matrices). In a world of massless quarks, the 
cost of producing a quark-antiquark pair with total angular momentum and momentum zero, 
is small. The vacuum can be seen as containing a condensate of qq pairs with strong attractive 
interactions. Mixing of chirality ()3U() implies that up and down quarks can acquire an effective 
(constituent) mass by moving through and interacting with the vacuum (see for example jZj). 

We can evaluate current matrix elements by differentiating the classical action 5*2 = 
/ d'^xC2 with respect to external sources. This gives the lowest order result: 

(0| J7'^75n|7r+(p)) = (0| |7r+ (p)) = iV2Fp'^ (32) 
(OI^VIO) = - (0| ^ |0) = -F^i^o (33) 

OS 

From the definition of the pion decay constant, 

iV2F^p^ = (0| dYl5U |vr+ (p)> 

we see that to lowest order we can make the identification F = F.,^. Relations ()32() & (|33() are 
only valid in the chiral limit, and are subject to corrections of order O (mq) . ()33() relates Bq 
to the vacuum expectation value. 

The combination Bouiq is experimentally determinable. The reason for identifying O (mq) 
O (p^) becomes clear if we expand ()29() to second order in meson fields. We then get the re- 
lations 

Ml+ = Bo{m^ + md) (34) 
M|+ = Bo{mu + m,) (35) 
M|o = Boima + ms) (36) 
Ml = i^o (m„ + + 4m,) (37) 

By elimination we obtain the Gell-Mann-Okubo consistency relation for the GS bosons 
SM^ = AMj. - l8|. This is reasonably^ weh satisfied using rriu — rud and M^g ~ 
(see 1] mixing). For higher mass resonances both linear and quadratic mass formulas give ac- 
ceptable relations [2]. This is because to first order in symmetry breaking we have 5 (m^) = 

^Putting in the numbers we find that we're oiT by approximately 20 MeV, but if we instead use hnear 
relations, the numbers get worse (~ 70 MeV). 
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(mo + 5m)^ — ttiq = 2niQ6m + ... But when expanding around a massless limit m and m? 
distinction becomes important. In the normahzed effective theory, the pion mass prediction 
is 

ml = {mu + TUd) Bo + {mu + nidf Cq + ... (38) 

However, there is no symmetry constraint to force the renormalized parameter Bq to zero, so 
the squared pion mass is (mostly) linear in the symmetry breaking quark masses. 

4.4 The O (p^) effective action 

Listing all possible operators invariant under discrete and continuous symmetries, transform- 
ing linearly, we can construct the tree-level effective chiral action of order p^. The number 
of terms can be reduced to a minimum using the equations of motion for O (p^) , partial 
integration, the unitarity of U and SU (n) n-depending trace identities. C4 for SU (3) was 
first determined by Gasser and Leutwyler |S] and is given by 

+L4 (d^UD^'U^) (xU^ + Ux^) + L^l^i^D^UD^'U^) (xU^ + Ux^)) 

+L6 [xU^ + Ux^y + Lr (xU^ - C/x^)' + Lg (xUhU^ + X^Ux^u) 

+iL9 (^Lf^^Dt^UD'^U^ + R^^Dt'U^D''U^ + (l^^UR^^'U^) (39) 

where Lj are the expansion coefficients that must be determined phenomenologically. Both 
C2 and £4 are invariant under 

U^U\ x^x\ ^i^^r^ (40) 

since the trace is invariant under cyclic shifting of matrices. This the so-called intrinsic 
parity operation (originally introduced by Witten [2j), which operates on the function, but 
not on the space-time coordinates. If = r^, then we can have terms containing an odd 
or even number of pseudoscalars, but no transition between the two, as this would violate 
intrinsic parity conservation. £^ is equal to if we have only EM interactions or direct meson 
interaction (no external fields). Including the W fields, the two sectors of differing intrinsic 
parity are coupled and intrinsic parity can be violated. Hence, both £2 and £4 are unable to 
describe anomalous processes like 3-k2K or Tr^'j'y. For further details see P]- 

What have we missed? There is the possibility of terms that transform non-trivially under 
G, but still preserve chiral symmetry if their variation under G is a total derivative .2 . . This 
type of term is anomalous. 

In order to describe anomalous processes we need terms that take the axial anomaly into 
account. The first terms that contribute are O (p^) and make up the so-called Wess-Zumino- 
Witten anomalous action. 

4.5 Anomalous Processes 

In this thesis we consider a number of anomalous processes. First we have the pseudoscalar 
to 77* decays: tt^/t] ^ 77*; where 7* is on-shell or off-shell going to an e^e~ pair. We also 
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Figure 1: Feynman diagrams for the anomalous processes. 

consider the semileptonic weak decays tt~^ /K^ — > 'ye'^Ve, vr^vr^e^i^, — > TT^TT^e^u 

and the three-pseudoscalar-photon interactions 77r~ tt^tt^ and ry — > 7r+7r~7. The r] cal- 
culations were performed for both the octet ryg and the singlet tjq component, but only the 
octet is relevant at the level we're working at. The Feynman diagrams for the processes are 
depicted in figure 1. 

5 T]' (958) 

Invoking the quark model, we see that the quantum numbers of ir^ ,'ir~ ,tt^ , , , 
and r]s are the same as for ud, du, {uu — dd),us, su, ds, sd and {uu + dd — 2ss). Extending this 
logic, we find one more state uu + dd + ss, identifiable with the next lightest pseudoscalar 
meson in turn, the r/' (958) . The rj' has an abnormally large mass compared to the other 
pseudoscalar s. This is because it receives a mass contribution from the axial U (1)^ anomaly 

jHj. By Noether's theorem the axial SU (3) singlet current is J^^J = Ylq=u d s'l^t^^^l^ 
divergence of which receives an anomalous contribution 

d>'4'J = '-tG%G'^^'' + 2im, q^,q ; (5"^^ = e^^^^G^^ (41) 

q=u,d,s 
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where G^y is the field strength tensor for the strong interaction and gi^'^'^l^ the antisymmetric 
Levi-Civita tensor. Taking the divergence of the vacuum to rjQ matrix element, we get 

{0\4'J\mip)) = iF,,p,e-^P-^ ^ (42) 
{0\d^4''J\7^o{p)) = F,,ml^ (43) 
K< = ir'^i^lG^G^'-nVoip)) (44) 

mq-^O ''0 " ^ 

which tells us that the rjQ mass is non- vanishing in the chiral limit. If it were not for the 
anomalous GG-term the U (1)^ symmetry would be approximately conserved and break dy- 
namically along with the chiral SU (3) , producing a nonet of GS bosons. 



6 Tjs — r]o mixing 

SU (3) breaking in the quark mass matrix causes the singlet rjQ and octet r/g components to 
mix, yielding the physical states rj and r/'. 



It?) \ _ / cos^ -sme\nT]8)\ ( Ivs) = ^ {uu + dd - 2ss) 
W) J ~ \ sine cose J \ \r]o) J ' | \r]o) = -L (^uu + dd + ss) 



(45) 



Quantum mechanically, other states with / = can also enter the mix. These we assume to 
be too heavy to be of significance. Another crucial assumption is that the mixing does not 
depend on the energy of the state, which allows for this simple phenomenologically motivated 
model |12j . Making use of the Tr^,rj,rj' to 77 data, the mixing angle is determined to be 
approximately —20°. With this angle we find the quark content 

T] ~ 0.58 {uu + dd) - 0.57ss (46) 
Tj' ~ 0.40 {uu + di) + 0.82SS (47) 

The ss is decreased for the r/, and increased for the r/', accounting for the large mass difference. 
Since the up quark has a greater charge magnitude, the r? — > 77 amplitude is boosted. 
Heuristically, this can be seen by 

sm6\ 9 / cose sin6'\ 9 / 2 cos 6* sin( 




= 

e ~ 20° 



(48) 



so that there is a factor 2 difference in width if we take the mixing into account. This mixing 
model can also be directly encoded into the chiral coefficients, but offers an explanation for 
lowest order experimental-theoretical discrepancies. 



7 Anomalies 

Anomalies are said to appear whenever a classical symmetry is violated by the existence of 
quantum corrections. They are crucial in expanding the theoretical framework, as they signal 
new physics in the standard model. In this section we examine the origin of the anomaly. 
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7.1 Classical vs. Quantum Symmetries 

Let's compare the classical Noether current with the one obtained from the full quantum 
theory using path integrals. A more complete treatment of the subject can be found in [H]. 
The infinitesimal transformation 

^p^ ^ = + £ (2;) f. ((^) ^ (49) 

where e (x) is the infinitesimal parameter (which is temporarily given a coordinate dependence 
for the purposes of Noether current construction) and fi is an arbitrary function of the fields, 
implies a Noether current and an invariance condition: 

dC 

{x) = ^ C {v\ dip') = C {ip, dip) + J^d^e (50) 

So that the Lagrangian is invariant if e is a constant. (jSUf) represents the classical symmetry. 

In the path integral formalism, all matrix elements can be obtained via the generating 
functional which is a functional of the sources. 

W \j] = e^^b'l = / [dip] e*/'^*^(^(¥''^v')-jV) (51) 



where the classical source field j (x) is used to probe the theory, and where [dip] stands 
for integration over all possible values of ip {x) at each space-time point. Matrix elements 
describing physical processes are obtained by taking the functional derivative^ of the logarithm 
of the path integral. 

(0| T {ip (x,) ...^ (x,)) |0) = (i)" ' ''A]^ (52) 

"3 \Xk) ■■■OJ \Xp) 

We can study the classical Noether current (x) by coupling it to a classical source field 
and inserting it in the generating functional (|51j) . 

W [Vf,] = J [dip] e' ^ d^^i^iv^d^)-^^Jn (53) 
(|52|) lets us take current matrix elements (denoted by bar), 

= ,54) 

dv/, [x) 

5\nW [Vf,] = InW [v^ + 6Vf,] -InW [v^] = -i J d^xJ^" (x) 5v^ (x) (55) 

If we choose = — (x) , then 

5e\nW[v^] = \uW[v^-d^e]-\nW[v^] (56) 
= i [ d^xJf" (x) df,e (x) = -i [ d^xe (x) d^Jf" (x) (57) 



■^Functional differentiation is defined by 
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where in the last step partial integration has been used. From this follows that if 5e In W = 
then all matrix elements are divergenceless, i.e. d^J^ {x) = 0. 

Since we are integrating over all values of ip {x) at each point in space-time, one can argue 
that it should make no difference if we shift the origin of integration at each point x by a 
constant, redefining ipi {x) = ip\ (x) — e (x) fi {ip), as in (|49|) . with an accompanying Jacobian 
J = l. 



h^W[v^-d^e] = 




= y [dv;^] e*/°'*^-(^(^''^^')-^'"^'') = lnT^[z;^] (59) 

This would imply d^J^ (x) = 0, in accordance with classical current conservation. However, 
as was first shown by Fujikawa shifting the origin like this is not always allowed. The 
transformation H49|) can have a Jacobian ^7 7^ 1, so that the path integral measure is changed. 
If the change of variables is non-trivial then d^J^ (x) 7^ 0, and we have encountered an 
anomaly. 



7.2 The [/(l)^ Axial Anomaly 

In the chiral limit, the QCD Lagrangian is invariant under the global U (1)^ axial transfor- 
mation Ip e~*^'^^V- Using Noether's theorem we get the classically conserved singlet current 
"^il! — Yl Qlfj-l^Q^ with d^jf^ = 0. However, in going through the full quantum theory anal- 
ysis, it is found that the divergence of the matrix elements of the axial vector singlet current 
is proportional an anomaly, 

d^jf^ = ^-^Gl,G^^'' ■G1,^e^''^PGlp (60) 

Important consequences of the anomaly involves a large anomalous contribution to the vr*^ — > 
77 decay rate, the prevention of r/' becoming a Goldstone boson (thus separating the chiral 
SU (3) spectra into one octet and one singlet part). 

There are two main approaches one can take in examining the anomaly: 1) a direct 
calculation via the U (1)^ — > gg triangle diagram, or 2) the path integral analysis. The direct 
approach was originally taken by the discoverers of the C/ (1)^ anomaly'^: Adler ^1], Bell 
and Jackiw ^3]. They calculated the matrix element for the [/(l)^ current going to two 
gluons via the Feynman diagram in figure 2. In the resulting four-momentum integral we can 
make a change of variables, yielding an expression that is compatible with conservation of 
the vector color current or the U (1)^ current, but not both. Phenomenologically, we know 
that only the vector color current is conserved. In the presence of the U (1)^ anomaly, the 
axial symmetry is not even approximately conserved. The anomaly is also present in the 
U (1)^ 77 diagram, producing equation (|Hn|) with as ^ oiem and G^jy — > F^y. 

The path integral approach was taken by Fujikawa ^3]) and clarifies the origin of the 
anomaly. In the path integral formalism we introduce the singlet axial current coupled to an 
axial current source into a functional of the gluon field A'^. 

W[a,,Al] = I [dV;][dV']e'^''"('''^^"('^''^'^^)""''-'^°'') (61) 

^A.k.a. the Adler-Bell-Jackiw anomaly. 
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Figure 2: U (1)^ — > gg diagrams leading to an axial vector anomaly. 



The crucial point is that in redefining the fermion fields to absorb the e term of eq. (|59jl. the 
path integral measure is changed: 



where the Jacobian J' is independent of the fermion fields. The Jacobian is divergent, but 
can be regularized (as was done by Fujikawa) by removing the high energy eigenmodes of the 
Dirac field in a gauge invariant way. The freedom to change integration variables in the direct 
approach, corresponds to the freedom in choice of regulator in the path integral approach. 

8 The Wess-Zumino-Witten (WZW) Anomaly Action 

In section 4.4 we saw that if we want to describe anomalous processes, we must construct an 
effective Lagrangian that violates intrinsic parity. This can be ensured by always including 
the totally antisymmetric Levi-Civita tensor ^t^'^'^l^^ the presence of which will preserve normal 
parity while at the same time violating intrinsic parity. The origin of this symmetry breaking 
lies in the axial anomaly. It affects photonic processes like vr" — > 77, but also hadronic ones 
like KK tt^it^tt~ . For reactions like these we need to construct a Lagrangian, which takes 
the axial anomaly into account. This section contains a brief sketch of how this can be done. 
The derivation follows the sigma model approach employed by Donoghue et al jH]. 

The history of the anomalous O (p^) WZW action follows a somewhat crooked path. The 
first effective action analysis of the anomaly was made by Wess and Zumino, who arrived 
at their expression by integrating the anomalous Ward identities pp. Their Lagrangian was 
given as a Taylor expansion, and was later given a geometric interpretation and rewritten on 
a compact form as a five-dimensional integral with four-dimensional space-time boundaries 
by Witten (his form did not conserve parity, but was later corrected). 

8.1 Anomalous Terms 

We emulate QCD behavior by introducing a color quantum number (with Nc colors), and we 
are using the n, d, s quarks so the number of fermions is set to three (all with mass M). All that 
is needed to derive the anomalous action are the correct symmetries. A convenient starting 
point is the sigma model, which incorporates the essential features of spontaneous symmetry 
breaking and chiral invariance. Making use of the representation independence theorem one 
can then proceed with the non-linear sigma model using exponential parametrization. This 
facilitates the extension of the formalism from 2 to 3 flavors. Dropping all terms irrelevant 




(62) 
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to the anomaly and imposing the unitary change of variables tpL C^i^L^ i^R C'^Rj one 
arrives at a La grangian for fermions of mass M coupled to axial and vector sources. 

C = {ip - M) ijj (63) 

where 

D^ = d^ + tV^ + iA^j5 ; Vf, = -^(^ed^^ + ^d^e^ (64) 

; \ = -i{C^d^^-^df,C^) (65) 

Our model contains no gluons, but according to the Adler-Bardeen theorem this will 
not modify the result. Keeping in mind the section on anomalies, we see that the change of 
variables induces a change in the path integral measure, and the Jacobian must enter into 
the effective action. 



The second exponent cannot produce s^'^"^ at O (p'^) , so the anomalous effect must be lodged 
in the Jacobian. Having the distinct goal of calculating the Jacobian in mind, we introduce 
a continuous parameter dependence via the transformation ^ — > = exp(zrA • (p/2Ftj) = 
exp {iTcp). Transforming infinitesimally in r induces a change 5 J in the Jacobian, 

V; ^ V' = + ^sr^R I m m = I wwdi^y-^'^ (6?) 

^ 5J = e 



dT 



= -2itr(^75 (68) 

r=0 



Due to the anomaly, the Jacobian is divergent and must be regularized. Making use of 
Fujikawas method of removing high energy eigenmodes in a gauge invariant way we take the 
limit 

tr(^75 = hmtr ((^756"^^-^-) ; = 9^ + iV^ + ^^^^75 (69) 



From the definition of vector and axial- vector currents (jHU with ^ ^ .^t) we have that 

p^'p^ = d^d^' + a ; d^ = 8^, + iVr^ + a^^Kl^ = d^ + T^^ 

■ c7 = -2i,^A^f + i[a^ + iK^,i^f]75 (70) 

We now make use of the heat kernel expansion from thermodynamics and take the limit 
e — > 0, to obtain an expression for the regulated anomalous action T {(p) . 

tr^75 = i I d^xTi (^-^Y^e^a^j J d^xTr {^-f,a2) => (71) 

r ((^) = -iln J + ... = ^ 1^ dT I d^xTr (^l^e^.^pA'^A'^A^A^,) + ... (72) 

where Nc comes from the color sum and the ellipses signify terms without correct symmetry 
{e^iyaf}) ■ The only way to integrate this expression in a closed form, is to Taylor expand each 
axial- vector current around r = 0, and then integrate to obtain a series of local Lagrangians. 
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What Witten did was to treat the r as a time-Uke fifth dimension X5, with r = 1 corresponding 
to normal space-time. By expansion one can prove that the final the result depends only on 
normal space-time. To describe processes other than direct meson interaction, we must include 
external electroweak gauge fields via and r^. This will alter the covariant derivative, which 
in turn affects the Jacobian. The result is a gauge invariant tree-level action describing the 
effect of the anomaly at O (p^) in the chiral expansion. 

The WZW action can be written on a form more suitable for calculations, by performing 
the r integration where possible and changing variables [B]. 

S[U,e,r]wzw = -^|da^^'^""(sfS,^S^EfS^) 

/ '^'^W {wiU,e,rr'''' -W{l,e,r)) (73) 

where 



WiU,£,r) 



(L ^ R) 



(74) 



where 



and 



U = exp 



.V2 

' F 



L R stands for 



yL . 



+ 



_2_ 

K- 



TT ' 



^,1 - 
yL 



+ 



V 



VE 



(75) 



(76) 



(77) 



8.2 Example of Calculation 

As an example of how to use (|73j) we calculate the amplitude for -k^ — > '^e^e~ . The five 
dimensional term does not contribute as it contains too many fields. The second term in the 
normal space-time integral is just there for mathematical consistency, which leaves us with 
W^yap- It's easy to see that the L ^ R operation leads to the same result (since £^ = r^). 
So we can calculate just the L terms and then multiply by 2. Below we have used partial 
integration and the antisymmetry of et^'^'^l^ . 
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2V2 



2 2 
= _i^6y2 (MQ2) df,A,daAf3 = -i2jd^A,daAfs 

^(^0^^^.) ^ i_L^e^-/3fe^e.fe^4 (78) 

where the tt^ — > 76^ e~ amphtude can easily be obtained by applying the Feynman rules for 
QED. £ is the photon polarization vector and k the photon 4-momentum. The tt^ — > 77 decay 
is historically important as it paved the way for the theory of the anomalous sector. This is 
where Steinberger in 1949 first observed the effects of the anomaly. 

In contrast with the effective Lagrangians, the only parameter appearing in the WZW 
action is Nc'^. This is because it is a prediction of the QCD anomaly structure. In accordance 
with Adler and Bardeen there are no radiative corrections. One can show that Nc is 
an integer by making use of the fact that the five-dimensional integral can only depend on 
normal four-dimensional space-time. 

Squaring the matrix element, summing over photon polarizations and inserting one half 
for identical particles gives 



\^WZW\ ~ Q 2 z?2 2 Z-^ ^ii'^u^a'^fB^lJ.'i^'a'lB''^ /i fc fi, ,< 

Pol 



{ e/e^ = (0,1,0,0) 
eVe = (0,0,1,0) 
I fc/A;^ = ^(l,0,0,±l) 

^ • • • = -2 X (2013) (2310) - 2 X (1023) (1320) ; (2013) = e2kQe\kl etc. 

Pol 

+ (2013)^ + (2310)^ + (1023)^ + (1320)^ = \ml =^ 

r,, = — ^^^ = ^44lX — ^^-^-7.73el- 
IGtt 167rm^ 97r^F^ 2 2 vr ^^^3 54^3^2 ^ 

In excellent agreement with the experimental value T^^ = 7.7 it 0.5 it 0.5eV |18j . This is an 
important test for the number of colors as well as the anomaly structure and chiral symmetries. 

In section 10 we follow through with the O (p^) contribution to vr'^ — > 77*, allowing us to 
solve for the chiral coefficients. 



9 Loops and Renormalization 

The anomalous sector is subject to meson loop corrections. QFT loop corrections to the Born 
amplitude lead to ultra-violet divergences in the form of polynomials in masses or external 
momenta. All non-analytical divergences must cancel. The UV divergences are removed 
by introducing a finite number of counterterms at a given order. Poles from one loops in 
dimensional regularization can be obtained by considering quantum fluctuations around the 
classical solution of the EOM. 

*One can also introduce the number of fermions Nf, which we here set to 3. 
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9.1 Power Counting 



In this section we consider the order at which a general Feynman diagram contributes. The 
distinction should be clear from section 4.2 where we are asking which tree-level terms will 
contribute at a given order. Here we follow a purely diagrammatic approach (as in |^) to 
arrive at Weinberg's power counting theorem ^5]. We are specifically interested in the order 
at which loop integrals contribute. 

Consider a general diagram with Ny vertices and Nn vertices from so that Ny = 
A^n- The dimensionality of the couplings is M^<^, where Nq = ~ ^) ^'^'^ ^ 

the characteristic mass scale. For Nj internal lines and Ne external lines we get M'^^'~^^ . 
There is a relation between internal lines, vertices and loops; 

Nj = NL + Ny-l = NL + Y,^n-l (79) 

n 

Remaining factors must composed of a power of energy times the logarithm of the dimen- 
sionless E'^//j? (where ^ is the renormalization scale). Putting all this together, we get the 
matrix element energy dimensionality 

^ - MN.+2N,+2^„N„-2 E''F {E/f,) ^D = 2 + Y,Nn{n-2) + 2NL (80) 

n 

So that a diagram with Ni^ loops contributes at E"^^^ higher than the tree-level used in 
the calculation. This simplifies calculations considerably, since at low energies, only a few 
loops need to be taken into account. Loop divergences are handled in the usual way, and 
can be removed by renormalizing the parameters of the theory. The general effective La- 
grangian, compatible with the symmetry conditions, must have enough parameters to absorb 
the divergences. 

The lowest order action for anomalous processes is already O (/) (the WZW action). At 
O [p^) we have the power counting: tree level diagrams from Swzw with one vertex, and the 
rest from £2- In this thesis we will proceed up to the O (p^) anomalous action. At O (p^) 
we have: tree level diagrams with one vertex from Swzw-, one from >C4 and the rest from C2-, 
tree level diagrams with one vertex from the anomalous O (p^) action and the rest from C2, 
one- loop diagrams with one vertex from Swzw and the rest from £2- 



9.2 Infinite parts 

Starting from a chiral-invariant effective Lagrangian, the divergences from the loop integrals 
are constrained by chiral symmetry. In the previous section we saw that one-loops contribute 
at O (p^) higher than the order of the Lagrangian from which they where calculated. The 
anomalous power counting laws tell us that chiral invariant counterterms in the anomalous 
Lagrangian of O (p^) will absorb the one-loop divergences, leaving behind the renormalized 
coefficients. To calculate the divergent parts of the £2 one-loops, we can expand around the 
classical solution to the EOMs of our chiral theory, i.e. the GS boson matrix. 



5S2 

dU 
u = 
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where is an Hermitian matrix that preserves unitarity. The one-loop divergences are ob- 
tained by expanding £2 (eq. (EHl)) to O (^^) : 

J d^xC2 = j d^xC2 + \ j d'^xCAije + O {f) (81) 

The divergent parts of the new term in (jHlj) can be pinpointed by calculating the second 
variation of and identifying the relevant operators. The power counting of section 9.1 
told us that the WZW action can contribute a vertex to one-loop diagrams. Including the 
WZW action, the equations of motion and the variation of second order in ^ now receives 
new contributions, and the corresponding anomalous quantities are calculated through the 
second variation of Swz- The calculation of divergent terms is very technical and the reader 
is referred to [3] for a more in-depth analysis. 

9.3 Meson One-Loop Corrections 

We begin by considering the PS 77* decays. If both photons are on-shell there are 
no infinite parts coming from the meson one-loops. This means that we should find no 
counterterms in the chiral coefficients describing the decay. If however, the photon goes to an 
e~^e~ pair, coefficients containing parts designated to cancel the infinities will appear. The 
meson one-loop corrections were calculated by Bijnens et al. [20] via the Feynman diagrams 
in figure 3. To fully describe situation we must include the results from the wavefunction 
renormalization and the decay constant corrections. For the semileptonic ir^/K^ ^e'^v 
decays, the one-loops corrections have been calculated in |2j. These have Feynman diagrams 
similar to those in figure 3. 

We are also considering the experimentally well charted 3-pseudoscalar-photon interac- 
tions rjTnr^ & tttttt^. The one-loop corrections were calculated in |2Z1 via the Feynman di- 
agrams in figure 4. Only on-shell photons are considered here. For Tnnr'y, k"^ is very small 
(<C m^) , and for r]7nTj A;^ = unless e~^e~ —>■ ^* —>■ 3it, which we do not consider. The decays 
— > 7r"'"7r~e"*~z^ and K" — > -K^-K^e^v are subject to similar loop corrections (as calculated 

In the loop corrections (listed in appendix A) we find the so-called chiral logarithms of the 
form vn^pg logrripg/ , where /i = rup is taken as the arbitrary renormalization scale. We also 
find the divergent terms proportional to A (see eq. H83() ). which will be exactly canceled by 
the counterterms at O (p^) . These are either proportional to pseudoscalar masses or particle 
four-momenta. The latter can be seen to vanish in the soft (low energy) limit. Most loop 
corrections were calculated using m„ = 7^ rus, and the Gell-Mann-Okubo relation can be 
used to eliminate the rj mass. 

10 Anomalous Lagrangian of O 

The construction of the anomalous O (p^) is a formidable task and was first completed for an 
arbitrary number of flavors {Nf) in 2001 by Bijnens et al. The only symmetries needed 
to construct the desired EQFT, are those of the initial and broken subgroup. One starts by 
listing all operators of a given order in the power expansion, compliant with all continuous 
and discrete symmetries. If we want to keep track of the phenomenologically relevant terms, 
it is imperative to reduce the list of operators to a minimum using all available constraints 
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Figure 3: One-loop diagrams for PS 77* 




and identities. The regularizing terms that cancel the infinite parts of the ultra-violet loop 
graphs contributing at O (p^) must also be calculated. 

10.1 Construction of the Effective Action 

The following building blocks can be used in constructing monomials of O (p^) in the chiral 
expansion: n^, /i^j^ = V^Uj/ + Vj/U^, f±fj.^, x± = u^X''^'^ ^^xu. The motivation for using these 
is that they appear in the calculation of the divergent parts. Also, the number of terms built 
of these can easily be reduced to a minimum. QCD symmetry conditions of parity, charge 
conjugation and hermicity must then be imposed on the list of operators. Since we are dealing 
with the anomalous sector the antisymmetric s'^'^^f^ enters, imposing further conditions. 

Use of partial integration, antisymmetry conditions, the mathematical Bianchi & Schouten 
identities and the EOMs reduces the number of monomials to 24 for the general Nf case. In 
the three flavor case we use SU (3) matrices obeying the Cayley-Hamilton relation, yielding 
an equation allowing for the removal of one more monomial. The Lagrangian density for 
A'^^ = 3 is then given by 

23 

Cf = Y,Cror (82) 

i=l 

where the monomials are listed in table 1. The infinities are removed using the MS 
dimensional regularization scheme subtracting both infinite and omnipresent terms: 

Cr = Cr + r?,^A ;A = ^-ln(47r)+7i.-l (83) 
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Z£'^-/3^X+ /-a/3]) 


5 
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£^'^"''(n^{V7/-7v,/-a/3}) 





12 




-6 









Table 1: O (p^) monomials & renormalization coefficients 

where Cf^^ are the renormalized coefficients, je is the Euler constant, /x the arbitrary renor- 
malization scale and d is the number of dimensions (— > 4). The coefficients rji can be deduced 
from table 1. Once the infinite parts have canceled, we can solve for the renormalized coeffi- 
cients Cj^^ by experimental comparison. The triple (anti)commutators in table 1 are defined 
as 

[a, b, c] = aba — cba 
{a, b, c} = aba + cba 



We put the non-abelian field strengths in the operator f!^" = uF^'^u^ ± u^F^" 



u. 



10.2 Infinite pcirts 

Calculating the second variation of the WZW & 5*2 action {O (^^) terms), and proceeding 
with dimensional regularization one obtains an expression for the divergent one-loop parts in 
arbitrary number of flavors {Nf) and colors {Nc) . 

Zi-loop - i67r2(d-4) 11527r2F2^^^i ^\ ^ ^ N] ) Nf' 

' b 6 ' 



+40f + + [-1 + ]^ j - - lOOlf - 30r4 + O 



+20Y, + OY, + ^Olf - 40S' - OZ + bOY, + ^OZ - -^OZ) (84) 
where the OY are the monomials listed in table 1. 
10.3 Example of Calculation 

We continue with the 7r° — > 76"*" e~ decay. The fastest way to work yourself through the list is 
by counting the minimum number of fields a monomial must contribute. Then many of them 
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can quickly be discarded. We are looking for terms that can accommodate one meson and 
two photons. By rewriting operators in terms of vector (v^) and axial vector (a^) sources, we 
can quickly see which ones will contribute (since only contains the EM field). The most 
useful tools in simplifying expressions is partial integration and the anti-symmetry of £t^'^°'l3. 

= ^^{mu- ma) e^'^'^^i^'^d^A^d^Af, 



O22 ^'^'^ ^23 contribute since the photon going to e+e is has fe^* / 0. However, 



7" 



to ncxt-to-lcading order in expanding the exponential, and F+^iy ~ 
so we can exclude O23. For O22 S^* 

O^^ = £'^^"^(«4V7>/+7-/+a/3}) 

= - ^e^^'^'-f^ {d^M{d^FL-y.FLap + Fl^pO^Fl^,)) 
= 16V2^e^'''"^ {MQ^)d^d^A,daAf3 
= ^pf^'^-t^TT'd^A^daAp 
Adding up the contributing monomials, we get the amplitude 

{8B0 (4m„ - ma) + 48Bo (m„ - ma) - 3k*^ C^) (85) 



The relations in section 4.3 can be used to rewrite the quark masses in terms of the pseu- 
doscalar masses. 

11 Models for the Chiral Coefficients 
11.1 Vector Meson Dominance (VMD) 

VMD is a phenomenologically very successful model that has not yet been derived from the 
standard model. It is based on the fact that higher mass resonances always enter the theory 
through virtual effects. VMD makes a dramatic appearance in the pion form factor, where 
the Breit-Wigner shape of the p meson can be seen. 
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The exchange of higher mass resonances can be used to estimate the finite part of the 
chiral coefficients. Roughly speaking, VMD states that for most processes, the main dynamic 
effect below the chiral breaking scale comes from the exchange of a vector meson. Using the 
hidden symmetry formulation by Bando et al. 22;, and extending the formalism to include 
the anomalous sector J. Bijnens ^3, has constructed an anomalous vector meson Lagrangian. 
Vector meson exchange can be described by including the ideally mixed p-nonet 

P, = ^P.A^ + (86) 

in a chiral invariant way. Operator strings with the correct symmetry properties can be 
constructed using the covariant and hermitian building blocks p^y = df^p^ — d^p^ + ig[pf^, p^], 
^^^/if? ) where = expiM/F\/2. The coupling g can be determined from the /? — > tttt 

width. 

The heavy p-meson nonet, consisting of p^'^,K*^'^,K*^,uj & (f), will only enter virtually 
and can be integrated out, assuming that their masses are much greater than the momenta. 
The resulting expression will contain parameters that can be constrained by comparing with 
the radiative decay widths of vector mesons. We will not list the full anomalous vector 
Lagrangian here - the reader is referred to j3| for further details. As an example, coefficients 
can be chosen such that one obtains a Lagrangian describing direct ppM interaction (where 
M stands for a pseudoscalar meson). 

£ {p^ ^ p,M) = ^-|-eM-/3tr (d^p^d^pf^M) (87) 

This Lagrangian can be used to calculate the M 77* amplitudes via the pp resonance: 
M ^ pp ^ 77*. To obtain the connection with and one must first integrate out the p 
resonance. It should be noted that it is also possible to calculate the ppM amplitude using 
an ordinary Feynman diagram approach. Integrating out the vector mesons one then obtains 
the effective Lagrangian for M 77* processes: 

where is the EM field strength tensor, Q = ^diag{2, -1, —1) and S = exp{^^/2M/F). 
Similar O (p^) VMD expressions can be obtained for the other anomalous processes. 



11.2 Chiral Constituent Quark Model (CQM) 

The constituent quark model is based on the assumption that vertices has to come from 
constituent quark loops. This effect is merged with the chiral formalism by including the 
constituent quarks in a chiral invariant way. Here follows a brief introduction to the model, 
for a more in-depth and technical analysis the reader is referred to Bijnens [3j or Ball |23j . 
The euclidean space Lagrangian density can be written 

= Cqcd + Cm = qVq (88) 
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where 



V = -f^Df,+M ; D^ = d^ + igsG^, + ie^PL + ir^,PR (89) 
; 7W = -mg (c/tPi + t/Pij) 

where mg is the constituent quark mass. We neglect gluonic corrections and will consequently 
drop in (|89|) . This is a large distance QCD approximation that will reproduce the anomaly 
correctly and is expected to work reasonably well. A connection with the effective approach 
in terms of the GS bosons can be obtained by integrating out the quarks. 

^r{U/„r^) ^ J[dq][dq]e^'^'''^Q^^ ^ (90) 

r = logdetP = trlogP (91) 

and similarly for T* replacing V with . The determinant corresponds to a trace over 
color, flavor, Dirac indices and space-time. With T* and T we can obtain a real part and 
imaginary part T~ . It should be emphasized that this is only possible since we're working in 
euclidean space. In the physical Minkowski space-time imaginary parts are not allowed. Since 
we are ultimately interested in the anomalous sector, we evaluate the effect of the intrinsic 
parity operation, and find that — > ztF^. This means that T'^ contains an even and r~ 
an odd number of e'^'^'^l^. In order to preserve the ordinary and break the intrinsic parity 
symmetry we know that we must have an odd power of e^'^"^ . (|9H) gives 

r- = -logdetpl' (92) 

This can be manipulated to obtain the correct expression describing the anomalous sector. 
The process involves rewriting H92|) on a five-dimensional integral with four-dimensional space- 
time boundaries, singling out e'^''"^-terms and expanding in mq. All but the WZW terms will 
be suppressed by the constituent quark mass. In order to perform the integration over the 
extra time- like dimension (c.f. r in WZW section), one can then make a Seeley-DeWitt 
expansion, singling out the contributing terms with SDW coefficients Oj. 

In this thesis we want to test the CQM O (p^) predictions by comparing them to the 
chiral O (p^) coefficients fixed by experiment. We will thus focus on the O (p^) abnormal 
intrinsic parity effective action. The F"*" describing the normal parity sector, can also produce 
anomalous terms. This is because the EOMs have been used in rewriting F"^, and these 
contain both anomalous and non-anomalous parts. However, this type of terms can be shown 
to be of O (/) and are not of interest here. The final result involves the contribution of 
terms with Seely-DeWitt coefficients 03 (describing one external field and 3 PS mesons), 04 
(2 external vector fields and one PS meson) and as (5 PS interaction). The explicit form of 
Si=3 4 5 r~ (oi) is very lengthy and can be found in appendix B. 

12 Results 

12.1 Combining Amplitudes 

We are now ready to combine all the amplitudes: the WZW O (p'^) , the one-loops and the 
O (p^) terms. Following through with the example of the vr*^ — > 7e"'"e~ decay, we can verify 
that the divergent part of the loops are canceled by the parameters of O (p^) . This serves a 
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test for consistent calculations. Combining ((75|) & (|H5|) with the loop amplitude in appendix 
A we get 

[1 + ^.il^k'' - i^"(log 4 + log 4) + f 



512 
"9" 



y 

We see that the infinite parts cancel: = ■;^2i{ip-i |Afc*2 — ^^^ip-i iJ-'^~^Xk*'^ = 0. The function 
F comes from the evaluation of the loop integrals and can be found in appendix A. 



12.2 Theoretical Quantum Field Calculations 

The results of the theoretical calculations are displayed in table 2, where the first term in the 
brackets is due to the WZW action. For the sake of completeness the one-loop contributions 
All are listed in appendix A. 



Process: 


Amplitude: 


vr*^ — > 77 


[1 _ 128^2 (20m|, + ml- 15m2) C^^ 
-12if TT^ {Ami -ml- Sm^) Cf'] 


77 ^ 77 


[1 _ 128 (4^2^ _ 5^2 ^ 3^2) ^2cW + {ml - m2) vr^Cg^] 






rj — > 7e"'"e~ 


[1 _ 128^2 (12^2, - 15?7i2 + 9ml) CY' - 36if 7r2 {ml - ml,) C^^ 
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+ 4- 
TT^ — ^ 76^1^ 


[1 - ^vr^m^Cr + f vr^ (g^ + fe^) C^^l + A^^^,^. 


K''^ — 7e^z/ 


[1 _ 256^2^2^^-^^ + 2567r2 (m|, - ml) C^' + f tt^ (g^ + k^) C^r^ 


jTT^ — > 7r"'"7r~ 


[1 + GATT^ml + Cf"- - CYi'' - Cf^'') 


77 — ^ 77r"'"7r~ 


[1 + 1927r2 (m2 - m^) [C^^- - C^-) 
+f (4m2, - 3ml + 5"^') (^5^' + 2Cf 




+64^2^2 ^5^r _ 32^2^ ^ (^j^r _ i287r Vpk<^iT 




(p^o - p^o)[-167r2gCj^'- + 647rVCir'^ + 327r2p;^C}^n 



Table 2: Results of theoretical quantum field calculations. 



k stands for the photon and q for the di-lepton four-momentum. The other momenta have 



been labeled as needed. The Wf^ has been replaced with the leptonic current 

with Gp/V^ = g'^/8M^. In the M — > 77* decays, note the absence of loop corrections in the 
case of real photons. It is interesting to see that the kinetic part of the electroweak fields in the 
decays with one pseudoscalar and two external fields are all connected to the coefficient. 
From experiment we expect this chiral coefficient to be rather large with respect to the others 
in the decay amplitude. obviously becomes increasingly significant at higher energies, 

so it is important to extract a good value. In the cases where they are active, €^^,0^^^ 
and play a similar role - they are also connected to kinematical factors, but do not 

vanish in the soft limit. Inspecting the corresponding monomials in table 1, we can trace the 
kinematical dependence to the presence of an extra d-y 4-derivative. Other monomials have 
the same property but forcably contain a minimum of fields that exceeds the number allowed 
by the process. 

The Gell-Mann-Okubo relation has been used to eliminate the r/g mass in all processes 
except for those involving an rj. Terms connected to C^*" in 7r^je~^e~ and to Cj^*" in ir^je'^i' 
are proportional to m„ — m^, and do not appear in table 2 since we are working in the isospin 
limit. The loop corrections must be recalculated if we wish evaluate the amplitudes away 
from the isospin limit. 

In all processes containing an on-shell photon, the term proportional to has been kept, 
even though this will be set to zero when comparing with the experimental data. Such terms 
will become useful when we are comparing with the VMD &; CQM predictions. 



12.3 Experimental Comparison 

Experimental data in the form of slopes, form factors and decay rates will allow us to extract 
numerical values for the chiral coefficients. The slope parameter b is defined as 

i.e. the factor in front of the off-shell photon squared four-momentum, normalized by the 
on-shell amplitude. To calculate the width we make use of the standard formula for two body 
decays, except in the case of 77/71"'" vr" where it's necessary to perform three body phase space 
integration. The form factors allow for easy comparison, as they can be directly related to the 
matrix element. Below, all cases are treated individually - the results are summarized in table 
3. The coefficients have been solved for using the least square solver in the MINUIT program 
(part of the CERN programming resource library). The errors are mainly due to experimental 
uncertainty, but all error limits have a contribution coming from the decay constants, of 
which we still have relatively poor knowledge. Measurements on charged pion decays give 
i^TT = 92.4 lb 0.33MeV. Next-to-leading order values for the other decay constants can be 
extracted through wavefunction renormalization as in p]. Many sources use an alternative 
convention and quote /m = V^Fm- 



12.3.1 TT^/r] 77* 

For real photons we form the decay rate by squaring the matrix element and tagging it with the 
appropriate factors. The squaring involves inserting one half for identical photons and then 
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Process: 


Experimental input: 


S-oZwed Coefficients [lO-^MeV-^] 


vr — > 77 & 
7? ^ 77 


Width [IH] 
Width dH] 


CY'' ^ 0.013 ± 1.17 
Cg^^ ~ 0.76 ±0.18 


vr*^ 7e^e^ 


Slope parameter [2^ 


CY'' - 6.52 ± 0.78 


— > 7e"'"e~ 


Slope parameter [2^ 


CY'' - 5.07 ± 0.71 


vr^ — > 76'*' 


Form factor jl8j 


CY'' ^ 20.3 ± 18.7 


ii'"'" — > 76"'' 


Form factor ^H] 


CY"^ ~ -6.37 ± 4.54 


77r'^7r'*'7r~ & 


Form factor |2H] 
Form factors |22] 


CY'^ ~ -0.32 ± 10.4 
CY"^ ^ 0.28 ±9.19 
(jWr ^ 28.50 ± 28.83 
CY"^ ~ -74.09 ± 55.89 

C'Wr ^ r)Q QQ _|_ 1 1 1/1 

(jWr ^ _25.30 ± 23.93 




Width [IH] 


^Wr _ (jWr ^ 21.67 ± 17.41 


K+ 7r°7r°e+i/ 







Table 3: Solved coefficients 

summing over all polarizations as in section 8.2. Width data is quoted in ^H] as a weighted 
average of several measurements, yielding F^o^^ = 7. 836 ± 0.523 eV and r,,^-y = 465 ± 44 
eV respectively. Note that the calculation was performed with the Goldstone boson matrix 
(|11|) . i.e. using r/g. Instead of explicitly implementing the mixing model of section 6, we 
conveniently let the chiral coefficients supply the needed extra factor. Proceeding in this way 
we are left with two equations, allowing us to solve for CY^ and CY^- 

For on-shell photons we can extract using (jSHI). To this end, we use the results of the 
CLEO II detector differential cross section measurements [2H. Fitting the form factor data 
they arrive at the pole parameters A^o = 776 ± 38MeV (= 6"^/^) and = 774 ± 49MeV. 
The similar slopes indicate that the two meson wavefunctions are nearly identical. Looking 
at table 3 we see that the two extracted values for CY^ agree within the error limits. Also 
note that CY^ is an order of magnitude larger than CY^ and CY^- 
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12.3.2 7r+/K+^7e+i/ 

For these decays the Particle Data Group ^S] quotes constant vector form factors (Fy), as 
no momentum dependence can be seen in the region of the experiments. The axial current 
is not of interest here as is belongs to the non-anomalous sector. Experimental vector form 
factor measurements yield the matrix element 

M {SDv) = ^^^^e^'l" F^e^^^rk^q^ ; F{^ = (0.017 ± 0.008) (94) 
V 2mp 

Comparing with the pion amplitude in table 2 we find 

where the one-loop corrections from appendix A have been inserted. Setting the photon 
and dilepton four-momenta squared (/c^ and q"^ respectively) to zero, allows solving for Cj^^. 
The measurements are rather imprecise, so the extracted value for should simply be 

regarded as an upper limit. The M — > 77* value is closer to the truth. As for C^^^, this 
coefficient can only be extracted from the decay, or alternatively from the vr"*" decay 
recalculated away from the isospin limit. Considering the poor data presently available, the 
recalculation hardly seems worth the effort. Repeating the above procedure for the kaon, with 
F^ = (0.204 ± 0.070) , we can extract a value for Cy^^. The fact that this values is negative 
is not a problem, since there is nothing in the theory which contradicts negative values for 
the chiral coefficients. 

12.3.3 77r°7r+7r^ 

Experimental data [IHl does not reveal any kinematical dependence. The amplitude is express- 
ible in terms of the jStt coupling constant F^'^ = 12.9 it 1.4GeV~^. Parsing the amplitude 
from table 2 and the loop corrections from appendix A, we set = and the invariant 
Pij — {Pi ~^ Pj)'^ to their average value - one third into their respective kinematical range. 

-3.5ml < Poi = {PnO-Pn-f < ^ (pIi) ~ -Iml 
^fnl < Pm = (PttO - P^+f < 13m2 ^ {pl^) ~ 7ml 



j:pl = Emi^{ph) = -'i 



2 

m± 



This gives 



+4F {ml, -fml) + 4F {m^ 7ml) + 4F {ml, -Iml)} 
-GAiT^ml [C7|X" + C^'' - 2Cf " - C^']] 

producing one equation with the unknown combination {C^^^ + — 20^"^ — C^^) . What 
we have done is clarified in figure 5, where the theoretical form factor and F^"^ have been 
plotted as functions of pgi ^^'^ P12 ™ units of ml . We have simply aligned the form factor 
average with the experimentally observed, constant plane F^'^ . Note that the error on 
transcends the maximum deviation of the theoretical surface from the plane. 
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Figure 5: Experimental (dashed) & theoretical (solid) form factor. 

12.3.4 ?777r+7r~ 

To compare with available data in the form of the decay rate T^-yTTTr = 55.2 it 6.6 eV |18| . we 
must perform three-body phase space integration. We make use of the formula 

by rewriting the amplitude in terms of the invariants = ip+ and p^s = (P- + 

and evaluating the integral numerically. The final result is an equation with chiral coefficients 
4-6 and 13-14 as unknowns. This equation contains a term proportional to — . Since 
these coefficients only appear in r]'jTnr, they cannot be solved for individually, and all we can 
do is to extract a value for the difference between them. 

12.3.5 K+ vr+vr-e+i/ 

This decay provides the best experimental data in the form of an energy dependent vector 
form factor. Data from 4x10^ events where recently (2001) collected at the Brookhaven 
Alternating Gradient Synchrotron j25j over a broad kinematical range. The matrix element 
is quoted as 

M = ^V:,u (p.) (1 - 75) V (Pe) {V^ - A>^) (98) 

where we only concern ourselves with the hadronic vector contribution = He^'^P" L^PpQ^, 
with P = p+ -|- p„ , Q = p^ — p^, L = Pe+ Pu in units of rriK and where the dimensionless 
H \s a, function of the invariant dipion mass Mt^^^ = |p+ + p- \ . 6 datapoints for H are 
quoted in the range 280MeV < M^^-,^ < 380MeV. No angular dependence was detected in this 
energy range. Comparing ()98() with the amplitude in table 3, we can extract 6 equations with 
H(M^^, [pk — p+)^ , {pk ~ P-)^)- This requires rewriting the matrix element using relativistic 
kinematics and verifying that the dependence on (the polar angle of the vr"*" with respect 
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Figure 6: Experimental & chiral predictions of H. 

to the dipion in the kaon rf) is small. Let the z axis be parallel to the dipion flight direction. 
In the kaon restframe we have 

PK = {mK,0) 

where |p27r| — — M^^)^ /m^ with ~ 0. In the dipion restframe we have 

P2n = (^2^, 0,0,0) 

p+ = (i£;27r,0,p^sin6'7r,P7rCOs6'^) ;p, 
p- = ( i E2n , , -Pt, sin , -Pt, cos Ot,) 
PK = (-Ek,0,0, - IpkI) 

where we can extract Ek and \pk\ using the frame invariance of scalars: Px^'^pfn^ = 
Pk^^ p2n^ ■ ^^"^ recast the amplitude in the appropriate form and verify that the 
angular dependence is small. For example: 

PKP- = \ [mi+Mi,) - M_^yp^^;T;;;^cos0. (99) 

The second term becomes very small when Mtttt — rriK or Am-,^. In addition, the term is 
suppressed by M^yr and the cos factor. 

The 6 extracted equations, involving coefficients 2, 4, 5 and 13-15, are now merged with 
the equation from 737r, giving a total of 7 equations with 6 unknowns. This system is overde- 
termincd, which will help to reduce the errors. 

The data series has been plotted in figure 6, along with the prediction of the fitted chiral 
parameters. Evidently the data points at the high end of the spectrum deviates significantly 
from the prediction. It is unclear whether or not this is just a statistical fluke, or if it 
originates somewhere in the experimental setup. The dominating accidental background was 
from Ktt^'k^'k~ , with a pair detection along with an from the beam or coincident 

decay. However, this was reduced to 2.4 it 1.2% using a likelihood method. 

If we hypothesize that the deviating points are indeed due to a statistical fluke, this can 
be compensated for by extrapolating a line using the first four data points to the high energy 
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Process: 






C^"- ~ 0.78 ± 12.7 
Cf"- ~ 0.67 ± 10.9 
C^'' ~ 9.38 ± 152.2 
CiT' - -8-44 ± 69.9 
CJX'^ ~ 0.72 ± 15.3 
C^"^ ~ -3.10 ±28.6 


r] — > 77r^7r^ 


^H/r _ ^Wr ^ 4_g ^ 26.6 



Table 4: Coefficient values as a result of extrapolation, 
process: Amplitude: source: 

7r+ (i^+) ^ ^e+i/ ^(A;2 + g2) gj] 
^ttO ^ vr+vr- {p^^ + pg^ + pf^ + 3k^) ^ 



Table 5: VMD O (p^) predictions. 



region. Also, we exclude the 737r equation as the measurement is not very good. Proceeding 
in this way, we can refit the chiral coefficients, producing a radically altered slope (see figure 
6) and the CiS in table 4. The errors have been overestimated, reflecting the uncertainty in 
the extrapolation procedure. This is why the errors in table 4 have been inflated. 



12.4 VMD Comparison 

VMD is unable to predict the mass terms of the chiral theory, as the Lagrangian contains 
no explicit mass parameters. Table 5 shows the O (p^) VMD expressions that should be 
compared with the chiral predictions. The prefactors are the same as for the O (p^) chiral 
expressions, allowing for direct comparison with the terms proportional to kinematical factors. 

The M 77* contains no kinematical dependence for on-shell photons. For off-shell 
photons TT^ and rj give the same result: 

C^"^ = ^ - 8.01 X IQ-^MeV-^ 

2^ 64M27r2 

in excellent agreement with the result in table 3. 7r+ (K^) 'je'^v' offers no new constraints. 
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VMD 


ChPT 


ChPT (extrapolated) 


^22 


eill^ ^ 8.01 


J 6.52 ± 0.78 
\ 5.07 ±0.71 




riWr 
•-"13 


V2mj^^ 20.0 


-74.09 ± 55.89 


-8.44 ± 69.9 




2 128M^7r^ ~ "•'-^^ 


29.99 ± 11.14 


0.72 ± 15.3 


(jWr 


3 1 o nn 

2 128Af^7r2 - ^-LIU 


-25.30 ± 23.93 


-3.10 ±28.6 



Table 6: VMD & chiral predictions in MeV 
Two equations can be extracted from the ^3tt amplitude: 



From rj^TTTT we get: 



CV/ + Cr/ + CV,^ = (100) 



2Cr,^ - 4Cr/ + C^,^ = (102) 



p 

= mJp (i»'> 

Exactly the same two equations follow from Ke4. There are three independent equations 
since (|lU2p & (jlU3|) can be combined to give (|100|) . In fact, the situation is worse - (jlUip 
originates in mass terms produced by kinematical factors. In the derived VMD Lagrangian, 
approximations have been made removing some of these terms. This means that (|101|) is 
subject to corrections that will alter the numerical value somewhat. Only H102|) &: H103|) are 
exact predictions of the VMD model. 

S/GAM^TT^ PS 8.0 X 10-9 MeV-2 and evaluating the lefthand side of (ITO & (ITO using 
the values in table 3, we obtain 2C]f - 4C}X'' ± Cjf^ ^ (-244. 7 ± 148. 4) x lO^^ MeV'^ and 
2C|X'' - C^/ « (134. 1 ± 78. 17) x lO"'^ MeV'^. Using the results of extrapolation in table 
4 gives (-17.52 ± 188.3) xlO^^ MeV'^ and (9. 88 ±100.5) xlO'^ MeV'^ respectively - in 
much better agreement with the VMD prediction. 

Including (jlOlj) . we can solve for coefficients 13-15. The coefficient values that can be 
extracted using VMD comparison and the respective chiral predictions are displayed in table 
6 in units of 10"^ MeV'^. 

The values for C^"^, C^^^ & Cj^^ are not great, but we have to remember that coefficients 
13-14 also has a part proportional to mass terms that VMD cannot predict. By contrast, 



vr 
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is only proportional to k'^ and there VMD does fine. The values obtained by extrapolation 
are compatible with the VMD predictions. 

Of course, the larger the momenta is in the experiment, the more dominant the terms 
with coefficients connected to the kinematical dependence will be, making the VMD parts 
increasingly significant. 

The shortcomings of the VMD model is illustrated by the plots in figure 7, showing the 
respective form factors normalized to 1 (with the exception of Ke4,). The top left plot is for 
the TT'^'je^i' decay, where the upper plane represents the VMD prediction and the lower one 
the result from ChPT. Here the infiuence from the mass term proportional to is very 
small and C22 dominates, so the VMD prediction is well within the error limits. 

In the top right plot (for the K'^'je'^u), the effect of mass terms is more dramatic. Both 
coefficients Cj^^ & C^"^ are active and bring down the ChPT prediction considerably. The 
middle plane shows the effect of removing these terms, aligning the ChPT with the VMD 
prediction. 

In the lower left graph (showing vr*^ (r/) 77*) we see that the effect of mass terms is hardly 
noticeable for the pion decay. Removing the Cj^^ term from the 7r^77* amplitude does not 
produce a visibly different result. In the 7777* decay, both Cj^^ and Cg^"^ are active, and 
become augmented by the eta mass. Removing these shifts the curve considerably towards 
the VMD prediction. The error bars represent the error due to the kinematical term. From 
looking at the graphs, we can conclude that the pion processes -K^jeu &: 7r^77* are well 
predicted by VMD, as they are relatively unaffected by the mass terms. Similarly for j3tt &: 
77727r, VMD does better when only pions are involved. This is because the presence of the r] 
(or for that matter) induces a greater mass term contribution. 

The lower right graph is for the Ke4 process, where the chiral prediction (solid line) has 
been extracted using equations from i^e4 & ^3tt. Ke4^ produces 6 equations and is the major 
contributor to the predicted form factor. Removing mass terms shifts the ChPT towards the 
VMD prediction. The solutions in table 3, using all available data, does not produce the same 
slope as VMD. The possible statistical fiuke in the high end of the spectrum is responsible for 
this, and can only be eliminated by future experiments. Ignoring the fiuke by extrapolation, 
produces a slope similar to VMDs. And if we then remove the mass terms, we get almost 
perfect correspondence. 



12.5 CQM Comparison 

With the Lagrangians of appendix B, the corresponding O (p^) results from the chiral con- 
stituent quark model can be calculated and subsequently compared with the ChPT predictions 
of table 3. First of all, we need to extract a value for the constituent quark mass - mg. We 
proceed as in |3l, making use of the M — > 77 prediction and experimental slope parameters 
pij . Setting Rfj, = = eQA^ in T~ (03) , then performing the calculation and adding lowest 
order (WZW) amplitude gives: 
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Figure 7: Form factors. 



where C^^ = 1 Sz C^g = 1/ V^- Taking an average of the pole parameters from [25 and using 
to calculate the slope, we get 



A^o = 776 ± 38 MeV 
= 774 ±39 MeV 



(A) = 775 ±39 MeV 
mg = ~ 224 ± 12 MeV 



12m^ (A)2 Vl2 

Comparing with the O (p^) expression in table 5 gives tuq = rrip/ \/T2 ~ 222 MeV, in excellent 
agreement with the empirically derived value. Since we are using the same experimental data 
to fix mq that we use to predict the chiral coefficients, we cannot compare the two approaches 
in the M — > 77 case. Using thq = mpj^fvi^ will of course give the same 0^2^ prediction as 
VMD. Note that no mass parameters appear in ()1U4() . Unfortunately, in deriving the CQM 
Lagrangian, some mass terms have been approximated away, as they do not alter the end 
result notably. CQM is in principle able give us predictions for all mass terms. We just 
have to remember that the equations derived by comparing with these are subject to small 
corrections. 

Table 7 shows the CQM O (/) amplitudes, with the same prefactors as in the ChPT case. 
The TT^ and the decays give the same (exact) prediction for C^"^ = l/512mQ7r^ ~ 
3.94 X lO-^MeV-^ It seems that the VMD prediction overestimates the C22^ value, whereas 



35 



Process: 


Amplitude: 

12m^ 


Lagrangian: 


ir^ (t]) — > 77* 


r- (as) 


TT"*" — > 76"*" 




r- (as) 


X"'" — > 76"*" 1/ 


^(-m^++3(z2 + 3/c2) 


r- (as) 


vr'^ ^ 77r+7r~ 




r~ (04) 


r/ 77r+7r~ 




r~ (04) 




3Qi^[m2^ ml 6qp++3q^] 


r~ (a4) 



Table 7: CQM O (/) amplitudes. 



CQM underestimates it by roughly the same amount. Both models do equally well with 
respect to ChPT. 

Comparing the mass terms we get the (inexact) value C^'' ~ 5.10 x 10^l°MeV~^ from 
the TT"*" decay. This can then be used in the amplitude to get C^'' ~ -1.44x lO^^^Mgy-^^ 
which is a bit to small to be taken seriously. 

Comparing kinematical terms, ij^tttt & Ke4 both give the same two equations. 



2C^/-C^/ = - \ ^ 6. 31 X 10-9 MeV-^ (105) 
5 liomgTT^ 

2C|^^ - 4CJX^ + C^"- = ~ 1. 26 X 10-s MeV-2 (106) 

5 128m|7r2 ^ ^ 

73-7r gives one equation that offers no additional constraints. Equations (|105|) & (|106|) are of 
the exact the same form as the VMD predictions (|102() &: H103() . and the predicted value of 
the RH side hes in between (8.0 x 10"^ MeV'^) those of CQM. 

Turning to the mass terms, we get one equation from 737r, which can be reconstructed 
with the two equations from ij^tttt. These are: 

-0C3>J- + 2CfI' + Cir = (107) 



with = — Cq^^'. Two additional equations follow from comparing mass terms in 

cr^cf'-c-:- = (109) 

If we wish to solve for the chiral coefficients using the (inexact) predictions from the mass terms 
(eq. (|107() - H110|) '). we must make additional assumptions. One approach is to parametrize the 
solution of the system of equations (|105jl - (|110|) in terms of (for example) & C]^^, and 
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Process: 


CQM 


CQM (extrapolation) 


tt"^ — > 76"*" 


C^'' ~ 0.51 ± 0.06 
C]^'" ~ 3.94 ± 0.43 




K'^ — > je'^v 


(jWr ^ _0.00143 ± 0.03 
Cs^'' ~ 3.94 ± 0.43 




737r, 7777r7r & 1^64 


C2^^ ~ 4.96 ± 9.70 
^3^^ ~ 5.07 ± 5.07 
~ 6.32 ± 6.09 
Cg^^ ~ 33.05 ± 28.66 
C}^'' ~ 14.15 ± 15.22 
CJX'' ~ 10.23 ± 7.56 
CJI''' ~ 19.70 ± 7.49 


C^'' ~ -0.074 ± 13.3 
C3^'' ~ -2.14 ± 6.54 
Cf''^ ~ -0.55 ± 9.05 
(^H^r ^ 34 51 _^ 41 ;l3 

C]!''' ~ -7.46 ± 19.62 
Ci^"^ ~ -0.58 ± 9.77 
C}!''^ ~ 8.89 ± 9.72 



Table 8: CQM & chiral predictions in MeV~'^. 

then vary it freely with the constraint tliat tlic solution should stay as close to the ChPT result 
as possible, given the error limits. Table 8 displays the results of applying this procedure on 
the full ChPT and extrapolated ChPT result, respectively. Looking at figure 8, showing the 
ChPT along with CQM and VMD predicted coefficients in units of 10~^ MeV~^, we see that 
the best result is obtained by extrapolation. The quoted CQM errors arise from the constraints 
in the fitting process and from the constituent quark mass error, and are not intrinsic to the 
model. 

Figure 9 shows the chiral and CQM form factor predictions for Tr~^^e^iy (top left), K'^je'^v 
(top right), 737r (bottom left) and Ke4, (bottom right). For vr"*' (K^) ^e'^u CQM does slightly 
better than VMD, but fails to predict the mass terms in the K'^^e'^v case. In the 737r process, 
CQM & VMD fit well with the coefficients obtained by extrapolating, less so with the full 
data set. In the K^i, both VMD &; CQM predict roughly the same form factor. 

13 Conclusions & Outlook 

To gain insight into the physics leading to the chiral Lagrangian, it is vital to chart the 
anomalous territory. Ultimately, we hope to arrive at a complete standard model derivation, 
capable of predicting all the chiral parameters. Going through the ChPT theory up to the 
anomalous O (p®) Lagrangian, we are faced with a large number of coefficients, that can 
be empirically determined through comparison with available experimental data. As much 
information as possible is extracted using the amplitudes calculated for a number of different 
anomalous processes. The extracted values then allow us to test the predictions of two other 
models for the C^^s : the vector meson dominance and the chiral constituent quark model. 
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Figure 8: Chiral coefficients (bottom: using extrapolation). 

VMD does very well in predicting the coefficients that are connected to kinematical factors, 
but fails to fully predict those that do not vanish in the low energy limit. Here the CQM does 
better, and gives reasonable predictions for all coefficients. 

To improve the accuracy in gauging the validity of these models, and to reduce the errors 
of the Cj^^s further, more detailed experiments must be made, probing the anomalous sector. 
Many such experiments are in the planning stages today, e.g. the Primex precision vr'^ lifetime 
measurements at Jefferson laboratories. 
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A One-loop Corrections 



Process: 


Amplitude: 


TT^ — > je~^e~ 


1 ^nval3^ u 1 

{-\k*\\og ^ + log ^) + f A;'^^ + 4 (^.2^ ^2) + p [k*\ml)]} 


rj 7e"'"e~ 


{-ifc*2(iog !!4 + log ^) + Wf.*2 ^ 4 (A;*2^ ^2) + F {k*\ml)]} 
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TT"'" — >■ 76+1/ 




— je~^u 


+4/ (P, ml, ml) + 2/ (g^, mj^, ml) + 21 {q^, mj^, m^) } 


77r°7r+7r'~ 


2 

(POI + P02 + P?2) log > + 1 (POI + P02 + Pf2) 

+4[F K,p2^) + F [mlpl^) + F (m2,p2^)] 


7?787r+7r~ 


3^{- (4m^ + Ipj,) log ^ + (4< - IpI,) log !j 
+|j.?2 + iF {ml,pl^) + If {ml,pl^)} 


iiT"'" ^ 7r+7r~e"'"i/ 


r 1 r 11 2i m2 ^ 2 i 2. 2-i rn^ 

+21 {{p+ +p-f, ml, mi) + 7 + p")' , m|^, m^^) 

+2/ ((p- +qf ,m\,ml] +1 ({p~ + ,mi,m^) 
+3/ (^^m2^,m2) +3/ (g^, m^^, m^)}] 




sin e ^uya0 1 ^ ^ ^/ 

[32;^{ / ((g+p)^m|,,rn2) 2^ ((<? + p)^ "^X' "^^) 
+/ ((g + pO' , rnl, mi) + \l ({q + p' f , m\, m^)] 



where 
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l.*2 

2 



I (fc^, m^, 777,2) = dx [mf — (ml — x — x{l — x)k^^^ x 



m? — (m? — mo) a; — x (1 — x) A;^ 
log ^ ^2 



(mf — ml — A;^)^ — 4A;^m2 



2 2 

,2^2 ^2 ^2 ^2 



iB {k'^,ml, ml) -\[m\ + ml) + ^k^} 



zi?(fc,mi,m2j - ^g^2(^l 2 log ^4 + 2P ^°Smi it^ u+ - J 



u± = y k"^ — {ml ± ml) 

B CQM Lagrangians 

5 pseudoscalars: 
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1 external vector field & 3 pseudoscalars: 

tr{2 (^D^r^^ + T,^D^i^^Y:^ D^^E^ Da^DpH^'E 

-3 (D^r^^ + s^D^^^^s) d^j:^d^j:d0j:^j: 

-2D^'EDaT,^ D^Dp'ET,^ + 2J:D^D^T,^ Da'BDpT,^ 

+3r,^„[D^DaT,^T,Dp'E^ D^'E + DaT,^ D^DpET,^ D^T, 

-Slut.lD^DaE'E^ DpED^T,^ + DaT,D^DpT,^'ED-yT,^ 

~2D^D^E{r,aE^ - St4„,SZ)^St} 
+2D^D^E^{i^aE - Sr^„,StZ)^S} 
+i2D^D^D^,'L[2D^D^Y)T,DpY) + 2Dfi'L^ED^Da'E^ 

~i2D^D^Di,Y)[2D^DaT,T) DpE + 20^^^'^ D^DaE 
+ED^DaT)DpE + DpED^DaT^Y]] 
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2 external vector fields & 1 pseudoscalar: 

tr{7{L>^r^^,r«^}L>,StS 

- [4/3, ^t] (pDjD^^'' + i^^^^^^^) } 
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